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Abstract 

We investigate the thermodynamics of a Fermi gas whose single-particle 
energy levels are given by the complex zeros of the Riemann zeta function. 
This is a model for a gas, and in particular for an atomic nucleus, with an 
underlying fully chaotic classical dynamics. The probability distributions of 
the quantum fluctuations of the grand potential and entropy of the gas are 
computed as a function of temperature and compared, with good agreement, 
with general predictions obtained from random matrix theory and periodic orbit 
theory (based on prime numbers). In each case the universal and non-universal 
regimes are identified. 
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1 Introduction 



Fermi gases confined to small volumes present distinct finite size effects in their ther- 
modynamic properties. Among the different corrections to the bulk behavior, some 
are smooth functions of the parameters of the system. They typically depend on 
geometric aspects of the volume occupied by the gas (e.g. surface and curvature 
terms). In contrast, other corrections are fluctuating functions. These quantum (or 
mesoscopic) fluctuations are related to the discreteness of the single-particle spec- 
trum, and are sensitive to the nature of the dynamics of the particles in the gas. 
In many cases they are small compared to the bulk, but may play nevertheless an 
important role (like, for example, in nuclear physics where shell effects are crucial in 
the determination of the shape of atomic nuclei [1, 2]). In other (more spectacular) 
situations, they provide the dominant contribution. This occurs for example in the 
electronic orbital magnetism in quantum dots - where they can be larger than the 
Landau susceptibility [3, 4, 5] ~, or in the persistent currents in mesoscopic rings 
[6, 7] ~ where the bulk contribution vanishes. In all cases the quantum fluctuations 
disappear when the temperature is raised. 

Instead of a detailed computation and description of the quantum fluctuations 
for a particular system, the aim here is to study their statistical properties. The 
interest of such an analysis is well known: using a minimum amount of information, a 
statistical approach allows to establish a classification scheme among the fluctuations 
of different physical systems. It also allows to distinguish the generic from the specific, 
and provides a powerful predictive tool in complex systems. A general description 
of the statistical properties of the quantum fluctuations of thermodynamic functions 
of integrable and chaotic ballistic Fermi gases, and of their temperature dependence, 
was developed in Ref . [8] . 

From a classical point of view, there are two extreme cases of single-particle 
motion, namely integrable and fully chaotic. The case of a mixed phase space (e.g., 
coexistence of regular and chaotic motion) being the generic one. In all cases, in a 
Fermi gas model the quantum thermodynamic fluctuations can be directly related to 
the fluctuations of the single-particle, discrete, energy levels. Schematically, there are 
two distinct types of single-particle fluctuations. On the one hand there are the local 
fluctuations, i.e. those occurring on scales of the order of the single-particle mean 
level spacing 6. For the two extreme dynamics mentioned above, these fluctuations 
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are known to be universal [9], namely Poisson for integrable motion, random matrix 
in full chaos. The second type of fluctuations generically present in any single- 
particle spectrum are long range correlations, that occur on a scale Ec, where 
is the energy associated with the shortest classical periodic orbit of the mean field 
(inverse time of fiight across the system). In contrast to the local fiuctuations, the long 
range modulations depend on the short-time specific properties of the dynamics and 
are therefore non-universal, either for integrable, mixed, or fully chaotic dynamics 
(thought their importance varies in each case). 

Properties related to both types of single-particle fiuctuations have been iden- 
tified in the nuclear behavior. For example, shell effects in the nuclear masses are 
a clear manifestation of the long range correlations, whereas nearest-neighbor spac- 
ing statistics and the distribution of widths of neutron resonances illustrate universal 
local fiuctuations. However, a global and comprehensive picture of the quantum fiuc- 
tuations has not yet emerged. An important additional ingredient are certainly the 
variations in the nucleus of the nature of the dynamics as a function of the excitation 
energy. But even in the case of a scaling mean-field motion where the classical dy- 
namics is energy independent (like when putting a Fermi gas inside a hard wall box 
or billiard) the physical behavior of the gas can be quite complex and rich, and both 
types of fiuctuations (local and long range) can in fact manifest in different quantities 
related to the system, or in the same quantity at different temperatures [8]. 

One might think that a fully chaotic single-particle dynamics implies univer- 
sal quantum fluctuations of the thermodynamics of the gas, well described by some 
random matrix model (Wigner, two-body random, etc). This belief is not true in 
general. The reason is that the thermodynamic fluctuations are universal only when 
they are controlled by the local single-particle fluctuations. And this is not always 
the case. For example, it has been shown, in particular, that the statistical proper- 
ties of the ground-state energy of a many-body system are controlled by the short 
periodic orbits (at least in systems with a ground state well described by a mean- 
fleld approximation, like the atomic nucleus) [8]. The corresponding distribution is 
therefore non-universal. As a consequence, it makes no sense to model the ground- 
state fluctuations by some random matrix model, that possess no information on the 
speciflc mcan-fleld short time properties of the dynamics. 

Our purpose here is to illustrate the richness of the thermodynamic quantum 
fluctuations and to explicitly check some of the predictions by considering a particular 
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example of a Fermi gas with a fully chaotic classical dynamics. The system considered 
is a mathematical model. It is achieved by considering a spectrum of single-particle 
energy levels constructed from the imaginary part of the complex zeros of the Rie- 
mann zeta function. The Fermi gas is obtained by filling the single-particle energy 
levels with an average occupation number determined by the Fermi-Dirac distribu- 
tion. This model of a chaotic gas was introduced in Ref.[10], where general motiva- 
tions as well as dynamical analogies were given. If this Fermi-gas model is viewed as 
an element of the periodic table, following nuclear physics terminology we referred to 
it as the Riemannium [10]. Aside its mathematical interest, physically the Rieman- 
nium is an important model because it possesses all the typical dynamical properties 
and characteristics of a realistic chaotic gas, while it greatly facilitates the numerical 
and analytical computations (e.g., all the ingredients required in a semiclassical de- 
scription of its thermodynamic properties are known). This is quite unusual in fully 
chaotic systems, and therefore offers a rather unique model to test general results and 
ideas. 

The basic quantities to be computed are the probability distributions, at a 
given temperature, of the quantum fiuctuations of the grand potential and of the 
entropy of the gas. We consider a degenerate gas in the grand canonical ensemble; 
other ensembles (e.g., canonical) may be considered as well. In particular, it can be 
shown that to leading order in an expansion of large chemical potential, the fiuc- 
tuations of the energy of the gas at a fixed number of particles coincide with those 
of the grand potential. Our results therefore apply to the total energy. In a semi- 
classical approximation, the fiuctuations are described by the periodic orbits of the 
classical dynamics associated with the single-particle motion [11, 12, 2, 13]. In the 
Riemannium, the role of the periodic orbits is played by the prime numbers. The 
basic equations and the relevant energy scales of the Riemannium are introduced in 
§2. In §3 we analyze the fiuctuations of the grand potential. The corresponding dis- 
tribution is non-universal at all temperatures (including T = 0, the case considered 
in [10]). The asymptotic moments of the distribution, obtained for a large chemical 
potential, are well reproduced by some convergent sums dominated by the smaller 
prime numbers (i.e., short periodic orbits). Odd moments are non-zero, implying 
an asymmetric distribution at all temperatures. At a finite chemical potential, large 
prime numbers introduce corrections to the asymptotic moments. These corrections 
to the non-universal leading-order behavior are, in contrast, universal, and described 
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by random matrix theory. 

The entropy fluctuations are studied in §4. Their character strongly depends 
on temperature. For a chaotic dynamics without time reversal symmetry, which cor- 
responds to the Riemannium, these fluctuations present several remarkable features, 
predicted in [8]. For temperatures T <^ 6, the variance of the entropy fluctuations 
increases linearly. In a second regime, when 6 <^ T <^ Ec, the variance saturates to a 
constant value. This initial behavior of the variance (i.e. linear growth followed by a 
saturation), and more generally of the probability distribution of the entropy fluctua- 
tions, was predicted to be universal, namely identical for any chaotic system without 
time reversal symmetry (with a given single-particle mean level spacing at Fermi en- 
ergy). It is well described by the fluctuations obtained from a GUE random matrix 
single-particle spectrum. The shape of the probability distribution strongly depends 
on temperature, and converges towards a Gaussian as temperature increases. How- 
ever, this universal behavior disappears for temperatures of the order Ec or higher, 
where the statistical properties of the gas do not coincide any more with random ma- 
trix theory. The probability distribution of the fluctuations is now described by the 
short periodic orbits, which are system dependent. In this high temperature regime 
the typical size of the fluctuations shows an overall exponential decay. All these 
predictions are worked out explicitly for the Riemannium and compared with numer- 
ical simulations. We have found very good agreement between theory and numerical 
experiments in all the regimes. 

Although only two particular thermodynamic functions are discussed, they il- 
lustrate all the rich possibilities that may be encountered in other functions. Similar 
methods apply to them. 

2 Thermodynamic framework 

In the grand canonical ensemble the equilibrium properties of a Fermi gas are defined 
by the grand potential 



where T is the temperature (we set Boltzmann's constant ks = I), IJ' the chemical 
potential, and 



Q(//,T) ^-T jdE log[l + 6^^^-^)/^ P{E) , 



(2.1) 
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is the single-particle density of states. In the Riemannium, the single-particle energies 
are given by the imaginary part of the upper-half complex zeros s = 1/2 + iEj of the 
Riemann zeta function ({s), Ej = 14.1347, 21.0220, 25.0109, . . . (the values used in all 
our numerical simulations were obtained from [14]). They arc expressed in arbitrary 
units, which are also those of temperature (since = 1). From the grand potential, 
other thermodynamic functions can be computed. We are particularly interested in 
the behavior of the entropy of the gas, 

5(/.,T) = -(g)^ . (2.3) 

The number of complex zeros with imaginary part less than E is [15], 

n(£;) = l-^£;iog7r-^7mlogrQ-z|) -^/mlogcQ-^^) • (2.4) 

Expanding the T function and using the Eulcr product representation of ( in the last 
term, this expression splits naturally into a smooth plus an oscillatory part, as in the 
semiclassical approximation in quantum systems. The density of zeros p{E) = dn/dE 
is written p{E) = p{E) + p{E), with 

?W = -:^EE^<=os(£;rlogp). (2.6) 

p r=l P 

In the Riemannium, the first term corresponds therefore to the asymptotic average 
density of the single-particle energy levels. In the fluctuating part p, the sum is made 
over the prime numbers p. The comparison of this equation with the semiclassical 
Gutzwiller trace formula for the spectral density of a dynamical system [12], p{E) — 
Sp.o. ^p.o. cos(5'p.o.(-E')/^) (the sum is over the classical periodic orbits), shows that 
each prime number may be interpreted as the label of an unstable primitive periodic 
orbit p of a fully chaotic system without time reversal symmetry, of action Sp — ETp, 
period Tp — logp, Lyapounov stability Ap = 1, repetitions labeled by r, and h — 1 
[16, 17]. Because of this mapping, in the following we often refer to the prime numbers 
as periodic orbits. Putting aside the issue of the very existence of a classical system 
whose quantum eigenvalues coincide with the imaginary part of the Riemann zeros, 
it is the formal analogy with the semiclassical theory of dynamical systems that 
makes the Riemann zeros in general, and the Riemannium in particular, interesting 
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on a physical basis as a model to test general theories of quantum chaotic systems. 
Additional support for a dynamical interpretation comes from numerical simulations, 
that clearly indicate [14] that the sequence of the imaginary part of the Riemann 
zeros obeys the Bohigas-Giannoni-Schmit conjecture for chaotic systems without 
time-reversal symmetry [18, 9]: asymptotically, their statistics coincide with those of 
eigenvalues of the GUE ensemble of random matrices [19]. 

Replacing the decomposition (2.5)-(2.6) of the density of states into (2.1) leads 
to a corresponding decomposition of the grand potential. A similar decomposition for 
other thermodynamic functions is obtained by derivation. We are interested in the 
statistical properties of the fluctuating part of those functions, and in their temper- 
ature dependence. In the forthcoming analysis, the different relevant scales are the 
following. 

i) Largest scale: chemical potential // 

ii) Intermediate scale: shortest periodic orbit 

Energy : = h/r^i^ = 27r/log2 9.06472 , 

Time : r„u„ = log2 , (2.7) 

Temperature : % = E^I2'k^ ^ 0.459224 . 

iii) Smallest scale: single-particle mean level spacing 

Energy : 5 = 1/p = 2tx j log(/i/27r) , 

Time : = hjb = log(;u/27r) , (2.8) 
Temperature : Ts = S/2t:'^ = l/7rlog(/x/27r) . 

The intermediate scale is related to the shortest periodic orbit (or time of flight 
across the system at energy fi) which, according to the previous analogy, corresponds 
here to the period Tmin = log 2 associated with the prime number 2. deflnes the 
outer energy scale of oscillations in the density of states (cf Eq.(2.6)). Due to the 
independence of the periods with energy in the Riemann dynamics, this energy scale 
is a fixed number. In the definition of the mean level spacing 6, we have used the 
asymptotic leading order approximation of the average density of states, is the 
so-called Heisenberg time. The 2tt^ factor (approximately 20) that relates energies 
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and temperatures is dictated by semiclassical arguments [8], and it has to be taken 
into account in numerical comparisons. A final important parameter is the ratio of 
Ec to S, which measures the number of single-particle states contained in a scale E^, 



g ^EJ5^ log{tx/27T)/log2 . 

In the semiclassical limit // — > oo, the three energy scales are well separated, S <^ 
Ec <^ II, and g oo. 

3 The grand potential 

Replacing p{E) by p in Eq.(2.1), integrating by parts, and applying Sommerfeld's 
approximation (valid when T <^ //), we obtain the usual thermodynamic relation 

n(Ai,r)=no(A^)-y p(/x)r2, (3.1) 

where Qq is the grand potential at zero temperature. For the Riemannium 

Uo{E) = log f _ 7 log^ , 

°^ ^ 47r ^ \2t:J St: 8 487r ^ ^ 

Terms of order up to E~^ in p have been included in the integration. The determina- 
tion of the constant c is of numerical importance when computing the fluctuations. 
Prom the work of Selberg on the function Si{t) [20] (which, up to an overall sign, 
coincides with the fluctuating part of the grand potential at zero temperature), we 
can extract its value 

^ ( fE r2 / TT- \ Q 7 ^c,(rE^ 

c=- log|C(a;)|(ix-lim / n(x) dx - ^log ( —] + —E^ - -E - , 
7rA/2 ^1^^ ^1 E-.oo\Jo ^' 47r ^V27ry 87r 8 487r J ' 

where n{x) is the exact smooth part of the counting function, computed from equation 
(2.4). Evaluating numerically the integrals, the value of the constant is found to be 
c ^ 0.75575. 

Replacing p{E) by (2.6) in Eq.(2.1), the oscillating part of the grand potential 

is 

ii(M.r) = -iE£^^f^^-(MrlogP). (3.2) 

The function 

k{x) — xj sinhx 
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takes into account temperature effects [2], and cuts exponentially the contribution of 
large prime numbers. If T > Tc, the contribution of each orbit, including the shortest 
one, is exponentially small. 

Equation (3.2) is an oscillatory function of fi that describes the quantum fluc- 
tuations of the grand potential. To study their statistical properties we define the 
average 

for any thermodynamic function /(//) associated to the Riemannium. The size of 
the average window A// has to be much smaller than // (in order to extract local 
properties), but large enough to contain a sufficient number of oscillations to make 
the statistical approach meaningful, i.e. <^ A// ^ 



3.1 Variance 

The most basic feature of the probability distribution of the fluctuations is the vari- 
ance. The average of the square of Q, can be written as an integral [8] 

mi^,T) ^ j^dr K{r) . (3.4) 

The function 

K{t) = ( 1] A^j cos[//(rj logpi - Tj \ogpj)\5{T - ^) ) > (3-5) 

is the form factor, i.e. the Fourier transform of the two-point correlation function of 
the single-particle spectrum; Ai — logpi/p["^^, and f = (r^. -\-Tp.)/2. This function is 
central in our analysis. In the semiclassical limit ^ oo it is conjectured to coincide 
with the GUE form factor 

KauEir) = min{T,Tjf} . (3.6) 

If we naively replace in (3.4) K{t) by Kgue{t), the integral diverges at r — > 0. This 
is because at finite g random matrix theory provides a wrong description of the short- 
time behavior of the form factor. For low values of r the diagonal terms in Eq.(3.5) 
should be used [21], 

Kn{r) = j:A^5iT-n) . (3.7) 
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This shows that the form factor is non-universal at short times, and that K{r) = 
for r < Tmin. The form factor behaves hke Kque only for times larger than some 
r*. The intermediate time r^, satisfies (Re7)^^ < t,, < Th-, where 7 is the classical 
resonance (eigenvalue of the Perron-Frobenius operator) with the smallest real part 
[22]. 

Since short times dominate the fiuctuations of the grand potential, in the limit 
r — s> we use in (3.4) the diagonal approximation of the form factor, and neglect for 
the moment the contributions of longer times. 

This is a convergent sum for all temperatures, and is independent of the chemical 
potential. At T = it gives the value ^ 0.079290 [10]. Since the temper- 

ature factor is a monotonically decreasing function of r, the variance exponentially 
decreases with increasing temperature, and the role of the short orbits (small primes) 
in the fluctuations is further enhanced. Fig.l compares Eq.(3.8) evaluated at different 
temperatures to numerical results obtained for the Riemannium (for reference, in the 
chemical potential window used to compute the fluctuations we have Ts 0.013). 

The agreement with the leading order description of the variance is excellent. 
However, as the chemical potential is lowered deviations from the asymptotic behavior 
are observed. The corrections may be obtained from Eq.(3.4), (il^) = {^'^)b + (f^^)off, 
by using the GUE form factor (3.6) for times r > r*. This gives. 



^ 2^ L ■ ^^-^^ 



As temperature goes to zero, the temperature factor tends to 1, and the integral is 
easily performed giving {Q^)oS — — l/(127r^T|). This provides a finite-// correction to 
the variance, which decreases as (log//)~^. For increasing temperatures the correction 
is smaller. The off-diagonal terms thus provide a universal correction to the leading 
non universal diagonal term (namely, the only specific information on the system 
that enters the correction is the average density of states at Fermi energy, through 
T„ = hp). 
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3.2 Higher moments and distribution 

We already showed that to calculate the variance the diagonal approximation is ac- 
curate because the short orbits dominate the fluctuations. A generalization of the 
diagonal approximation allows to evaluate all the moments of the probability distri- 
bution of the grand potential by a method developed in [8] and [10]. The mechanism 
is an interference process between repetitions of primitive periodic orbits. 
Defining the amplitude 

1 kItt T r log p) 

AAT) = -TT- \ 1 , 3.10 

27r H p^'^ logp 

the third and fourth moments of the distribution of VL are found to be, in the diagonal 
approximation, 

oo 

{Vl ) = 6 >^p,ri >^p,r2 >^p,ri+r2 ) (^-H) 

P ri,r2=l 

oo oo 

{VL ) = 2 ^ 1^ 4 ^ ■Ap,ri •^p,r2 •^p,r3 •^p,ri+r2+r3 ~ ^ ^ ■^p,ri ■^p,r2 
P ri,r2,r3=l ri,r2=l 

oo ri+r-z — l 

-\- 3 ^ ] ^ ] >^p,ri >^p,r2 '^p,r3 '^p,ri+r2—r3 ~l~ 3 (Q, ) . (3.12) 

ri,r2=l r3=l 

All the sums are convergent. We see explicitly that the third moment and the excess 
of the fourth moment (the term 3 (fi^)^ being the fourth moment of a Gaussian 
distribution) are different from zero. The distribution of the fluctuations of the grand 
potential is therefore a non universal asymmetric function that strongly depends on 
the small prime numbers. It is displayed, for different temperatures, in Fig. 2 [23]. 

Higher moments A; > 5 may also be computed, with increasing complexity of 
the result. Corrections to the diagonal terms coming from off-diagonal contributions 
exist, but as already shown they are negligible in the limit — > oo. We thus ignore 
them. 

These results are confirmed to high accuracy by the numerical data from the 
Riemann zeros (cf Fig.l and Table 1). 
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4 The entropy 



Prom Eq.(2.3) the entropy of the gas is expressed as, 

Sin, T)^js (^) p{E) dE^Y.^ (^) , (4.1) 

where 

s{x) = log(l + e^) - — ^ . (4.2) 

1 + e ^ 

The function s{x) is locahzed around x = 0, and decays exponentially for |a;| ^ 1. At 
very low temperatures (T <^ T5), the function S{p,T) has peaks as a function of n 
centered at each single particle energy level (Riemann's zero) of width <^ 6. At these 
temperatures the entropy is therefore zero for fi ^ Ej, and takes the value log 2 when 
fi = Ej. As the temperature increases, the width increases (while the height remains 
constant) and the peaks start to overlap. At temperatures T ~ T5 and at a fixed /i, 
only the energy levels distant by a few mean spacings from the chemical potential 
contribute to the entropy. The fluctuations are governed by the local statistics of 
the eigenvalues, which are universal (e.g., GUE). In this regime universality in the 
entropy fluctuations is expected. In contrast, at higher temperatures T ~ Tc, when 
peaks separated from by a distance of order Ec start to contribute to the entropy, 
the universality will be lost because information on the scale of the shortest periodic 
orbit enters. We now substantiate this analysis by an explicit quantitative calculation. 

Deriving with respect to temperature Eqs. (3.1) and (3.2), the smooth and 
fluctuating part of the entropy are given, respectively, by 

S{^i,T) = yp(/x)T, (4.3) 
Si^^, = E E cos(/. r logp) . (4.4) 

p r=l ' ^ 

The function k'{x) — dn/dx decreases also exponentially for x ^ 1, producing a 
cut-off for prime numbers satisfying logp » 1/{'kT). When a; — > 0, k'{x) vanishes 
linearly. The maximum of the function is located at x ~ 1.6061. Prime numbers that 
satisfy rlogp ~ 1/{'kT) give therefore the main contribution to the oscillations. At 
low temperatures large primes (long orbits) are selected. Smaller primes (e.g., short 
orbits) contribute as the temperature is raised. 
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4.1 Variance 

Analogously to the grand potential, the variance of the entropy fluctuations may also 
be written in terms of the form factor, 

^S^^^lj-,r^^^Kir). (4.5) 

We briefly review here the general results obtained in Ref.[8], and check them with 
the Riemannium. 

The two different regimes of the form factor described in §3.1 spht the integral 
(4.5) into three different parts, 

(4.6) 

The sum corresponds to the non universal short-time regime, whereas the two in- 
tegrals correspond to the long-time random matrix behavior. These different terms 
dominate the integral at different temperatures. 

* Low temperatures: T <^Ts. In this regime the maximum of k' is centered at times 
much larger than the Heisenberg time Th- The dominant term in Eq.(4.6) is the last 
one. We can extend the integral down to zero with a negligible error. The variance 
of the entropy is 

(S^) ^ h p(l^) T^h T/S , (4.7) 

where 

/2 = vr^ r dx '^-^ ^ 1.51836 . (4.8) 

JO X'' 

In this initial regime the growth is linear, with a slope proportional to the density 
of states. No additional specific information on the Ricmann zeta is present in this 
formula. Eq.(4.7) reflects the discreteness of the single-particle spectrum, and is thus 
a very general result valid for arbitrary systems, independently of their dynamics 
[8]. The smooth part of the entropy grows also linearly with temperature, and is 
proportional to the density of states [24]. Compared to the mean value S, the typical 
size of the fluctuations \J (S*^) are large (of order (T/5)~^/^). In this regime the 
quantum fluctuations dominate. 

* Intermediate temperatures: Ts <ti T <^ Tc. Now the maximum of k,' is centered at 
times in between Tmin and Th- The dominant term in expression (4.6) is the second 
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one, where we have introduced the GUE diagonal approximation of the form factor 
Kd{t) = T. Neglecting the contributions of short and long times we extrapolate the 
integral from to oo, 



In this regime, the size of the entropy fluctuations are therefore insensitive to tem- 
perature variations, they saturate to a universal constant. 

* High temperatures: T ~ ^-nd higher. In this regime the diagonal approximation 
of the form factor is still accurate, but the short time (non universal) structure is now 
apparent. At this temperatures the entropy fluctuations are also sensitive to the fact 
that below Tmin = log 2 the form factor is strictly zero. The dominant term in Eq.(4.6) 
is the flrst sum. The sum can be extended to all prime numbers and repetitions with 
negligible error, and gives a variance 



For T ^ Tc, all the terms of this sum are exponentially small, and the fluctuations 
vanish accordingly. 

To have a global description of {S^) that interpolates between the different 
regimes described above a numerical evaluation of Eq.(4.6) is necessary. The result 
depends on fi through the Heisenberg time. 

We have checked these predictions by a direct comparison of Eq.(4.6) to a 
numerical computation of the entropy variance as a function of temperature. The 
result is displayed in Fig. 3 (for reference, ~ 0.013). There is an excellent agreement 
with theory for all temperatures. The initial linear growth and the saturation to 
a plateau arc amplified in the inset. The size of the fiuctuations almost reaches 
the theoretical prediction (4.9) at T ^ 2.5 Tg ^ 0.03. The expected intermediate 
plateau is however short, because the temperatures Tg and are not sufficiently well 
separated (due to the slow logarithmic decrease of 6, even at these large values of n 
we are not sufficiently asymptotic). The exponential decay is also well described. 

4.2 Distribution 

In the regime T <^Tc the previous results confirm that the behavior of the statistical 
properties of the entropy fluctuations are universal. They depend only on the struc- 




(4.9) 




(4.10) 
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ture of the GUE form factor and not on any specific property of Riemann's zeros. 
The universahty is not expected to be vaUd only for the variance, but more gener- 
ally for the full probability distribution. To check this, we compare the probability 
distribution of the entropy fiuctuations obtained from Eq.(4.1) using two different 
single-particle spectra {Ej}: 

a) the zeros of the Riemann zeta function, 

b) the eigenvalues of a GUE ensemble of random matrices. 

In the first case the probability distribution is computed by varying the chemical 
potential in a small window, whereas in the second by averaging over the Gaussian 
ensemble. Both probability distributions, computed at different temperatures, are 
plotted in Fig.4. At low temperatures both distributions are almost indistinguishable. 
Notice the strong sensitivity of the distribution to temperature variations. When 
temperatures of order Tc are reached (cf part d)), the universality is lost. For T ~ Tc 
the fiuctuations are dominated by short orbits, and are system specific. The moments 
of the distribution can be computed by the same techniques used in §3.2 for the grand 
potential replacing Ap,r in Eq.(3.10) by Ap,r{T) = k^tt T rlogp)/(2r p**/^). 

5 Concluding remarks 

The use of the Riemannium as a test model in quantum mechanics is justified by 
two main reasons. First, by its genericity: the Riemann spectrum possesses all the 
generic features of a classically chaotic quantum system with no time reversal sym- 
metry. Second, by its practical advantages, namely all the necessary quantum and 
semiclassical information required to work out accurate computations and compar- 
isons is available. Hence, though this "number theoretic" model may seem somewhat 
remote from a realistic system, it provides an excellent arena to verify the non trivial 
quantum mechanical properties of chaotic systems. 

Based on semiclassical techniques and random matrix theory, several aspects of 
the thermodynamics of a chaotic Fermi gas have been verified. An accurate descrip- 
tion of the probability distribution of the quantum fiuctuations of the grand potential 
(or energy) and of the entropy of the Riemannium were obtained. In particular, the 
universal linear growth of the entropy variance followed by a saturation, with a further 
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non universal exponential decay, were confirmed. The size of the saturation plateau, 
predicted in the regime <^ T <^ Tc, was relatively small (see the inset in Fig. 3). 
This is due to the slow logarithmic asymptotic convergence properties of the Rieman- 
nium (in spite of the large chemical potential used in the numerical simulations, we 
are not very deep in the semiclassical limit. In fact, for the window analyzed in the 
figures the number of particles is around 10^^, with g ^ 35. For an atomic nucleus or 
for electrons in a metallic grain, this value of g corresponds to approximately 40 ~ 50 
particles. In the Riemannium, in order to have Ts and Tc separated by a factor of, 
say, 100, n need to be of the order of lO^''). 

The high accuracy of the results obtained for all the quantities studied confirm 
the validity of the different approximations employed. The present theory therefore 
provides a solid ground to go beyond and test realistic systems. Of particular in- 
terest is the interplay between mean-field approximations, residual interactions and 
dynamics. Some encouraging results in this direction were already obtained in the 
study of nuclear masses [25]. 

This work has been supported by the European Commission under the Research 
Training Network MAQC (HPRN-CT-2000-00103) of the IHP Programme. 
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Figure 1: Variance of the grand potential fluctuations of the Riemannium as a function 
of temperature. Dots: numerical results computed in the chemical potential window jj. = 
(267653402147 ± 6000) containing the zeros (10^^ ^ ^940) ^ (10^^ ^ 48684). Full line: 
theoretical prediction (3.8). 
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Figure 2: Normalized numerical histograms of the probability distribution of fl at different 
temperatures: (a) T = 0, (b) T = 0.1, (c) T = 0.3, (d) T = 0.5, (e) T = 0.8, (f) T = 1.2 
(same chemical potential window as in Fig.l). 
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Figure 3: Variance of the entropy fluctuations of the Riemannium as a function of temper- 
ature. Dots: numerical results computed in the same chemical potential window as Fig.l. 
Full line: theoretical prediction (4.6). Dashed line (inset): low temperature approximation 
(4.7). Dot-dashed line (inset): saturation value (4.9). 
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Figure 4: Normalized probability distribution of S at different temperatures. Full line: 
numerical results for the Riemannium computed in the same chemical potential window as 
Fig.l. Dash line: numerical results from a GUE single-particle spectrum, (a) T = ~ 
0.013, (b) T = 4 0.052, (c) T = 0.1, (d) T = 0.5, (e) T = 0.8, (f) T = 1.2. 
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T = 


T = 0.3 


T = 0.5 




(7.928 ±0.002) X 10-^ 
7.9290 X 10-2 


(5.9885 ± 0.0025) x 10-^ 
5.9886 X 10-2 


(4.2953 ±0.0015) x 10-2 
4.2953 X 10-2 




-(5.78 ±0.02) X 10-3 
-5.7822 X 10-3 


-(3.44 ±0.02) X 10-3 
-3.4377 X 10-3 


-(1.65 ±0.01) X 10-3 
-1.6508 X 10-3 




(1.480 ±0.002) X 10-2 
1.4814 X 10-2 


(7.625 ±0.005) X 10-3 
7.6273 X 10-3 


(3.586 ±0.003) X 10-3 
3.5869 X 10-3 



Table 1: Moments of J7 at different temperatures. The upper values (with the errors) 
were obtained from the numerical distributions of figure (2. a), (2.c) and (2.d). The lower 
values are the semiclassical results (3.8), (3.11) and (3.12), respectively. 
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